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FOREWORD 


This document presents results of a recently completed joint Boeing- NASA program to 
study the effects of winglets on flutter characteristics of twin — engine transport type 
wings and to verify flutter analysis methodology. This document is one of the two 
proposed NASA publications dealing with this study and contains details sufficient to 
permit independent vibration and flutter analysis. A second publication, a NASA 
Technical Paper (TP), is planned for 1985, and will contain a technical summeuy. The 
present document is in two volumes: 

Volume I - Low-Speed Investigations 

Volume II - Transonic & Density Effect Investigations. 

The two volumes are arranged such that each volume may be used independently of the 
other volume. The foreword and introduction are common to both volumes and are 
Included in each volume along with a complete table of contents covering both volumes. 


Mr. C. L. Ruhlin of Configuration Aeroelasticity Branch of NASA Langley Research 
Center was the test engineer for flutter tests conducted in the NASA Langley 16' 
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1.0 Summary 


Flutter characteristics of a cantilevered hi^ aspect ratio wing with winglet were 
investigated. The configuration represented a current technology, twin- engine airplane. 
A low-speed and a high-speed model were used to evaluate compressibility effects 
through transonic Mach numbers and a wide range of mass-density ratios. The results of 
the investigation are described in two volumes of this NASA CR and summarized in a 
forthcoming NASA TP. The results from the test in NASA Langley 16' Transonic Dynamic 
Tunnel (TDT) and analysis- test correlation are included in tills Volume II. 

The low-speed model was retested ir.- TDT to determine altitude or mass- density ratio 
effects. This model had been earlier tested In General Dynamics, Convair Division, San 
Diego low-speed tunnel, and the results are discussed in Volume I. The mass-density ratio 
was varied in TDT by testing the low-speed model in both air and freon. The 
configurations with winglet showed a switch in flutter mode, from nacelle vertical 
bending to second wing bending, due to decrease ui mass-density ratio. The mass-density 
ratio effects. Including the mode switch, were satisfactorily correlated between analysis 
and test. 

The high-speed model was tested in freon for a Mach range of about 0.6 to 0.91 and 
dynamic pressures up to 200 psf. Four flutter mechanisms were obtained in test, as well 
as analysis, from various combinations of configuration parameters. The winglet 
aerodynamic effects were significant and caused reduction in flutter dynamic pressure. 
The winglet related flutter for the configuration tested, was found to be amenable to the 
conventional flutter analysis techniques. The analysis showed that coupling between wing 
tip vertical and chordwise motions has significant effect under some conditions. 


1 




2.0 Introduction 


The interest in using wing- tip-mounted winglets to reduce drag for transport airplanes 
was stimulated by the work reported in Reference (1). One of the first applications of 
wlngHets was for the KC-13S airplane based on a potential drag reduction of about six 
percent estimated in Reference (2). The KC-13S Wlnglet Flight Research and 
Demonstration Program was formulated to design, fabricate and flight test a set of 
winglets to prove the drag reduction and other characteristics of the wlnglet concept. 
This program Included a low-speed wind-tunnel flutter model test and a flight flutter test 
program (Ref. 3). The critical mode during flight flutter test was a 3.0 Hz low- damped 
mode occurring with a light fuel loading at 21,500 feet altitude and with zero degree cant 
angle and -4 degrees incidence winglets. Plight testing for this configuration was 
terminated at 370 KEAS, rather than the test goal of 395 KEAS. due to low damping 
(g ^ 0.015). The low damping obtained for this mode was not predicted by flutter 
analysis. The lack of correlation was Judged to be due to limitations of current linearized 
aerodynamic theory and inability to represent transonic effects. Winglets have also been 
considered for the B-747 airplane as a part of the NASA Energy Efficient Transport 
Program (Ref. 4). Two flutter modes were obtained in the low- speed model test for the 
configuration with vdnglets. These flutter mechanisms were not present for the baseline 
configuration without winglets and were shown to result from wlnglet aerodynamics 
rather than mass effects. Flutter speeds for the configuration with winglets were 
significantly lower than the baseline configuration. It was suggested that the flutter 
mechanisms could be predicted by incorporating static-lift effects as with T-tail type 
flutter analysis. 

A transonic flutter model study of a supercritical wing with winglet for an 
executive-jet-transport airplane (Ref. 5) reported a good analysis-test correlation. The 
wlnglet addition decreased flutter speed by seven percent, of which a five percent 
decrease was due to the wing-tip mass effect. Thus, there was no significant reduction in 
flutter speed due to winglet aerodynamics. Results of another application of winglets for 
the DC-IO airplane, under the NASA Energy Efficient Transport Program, were recently 
published (Refs 6 and 7). A low-speed flutter model test showed that the winglets had 
generally detrimental effects on the flutter characteristics with small- to- moderate 



degradation in the basic wing flutter mode and a large degradation in a higher frequency 
wing flutter mode. During the flight test of the DC-IO airplane with winglets. SOO poimds 
of mass baluice was installed in each wing tip to ensure adequate flutter margins for 
flight testing. 

It appears from the available data that winglets generally caused degradation in flutter 
speed. The actual reduction in flutter speed varied with the configuration. The KC-135 
flight test experience of encountering an unexpected low-damped mode highlighted the 
technical risk Involved in flutter assessment of an airplane configuration with winglets. 
The only transonic wind-tunnel flutter test data available on a scaled airplane wing was 
for an executivejet-transport wing which showed a small reduction in flutter speed due to 
addition of a wlnglet. These considerations led to a joint Boeing/NASA program to 
develop a flutter methodology for winglet configured wings. A typical, current 
technology, twin-engine transport wing was selected as the basis for the study. A test 
program was outlined as follows: 

A. Pressure Model Test for Aerodynamic Data Base 

B. Low-Speed Test 

(i) Model Ground Vibration Test (GVT) 

(ii) Flutter Test and Parametric Studies 

(lii) Analysis-Test Correlation 

C. Test in NASA Langley 16' Transonic Dynamics Tunnel (TDT) 

(1) Retest of Low-Speed Flutter Model for Mass- Density Ratio Effects 
(ii) Selection of High-Speed Model Configurations 

(Ui) High-Speed Model GVT 

(iv) High- Speed Model Flutter Test 

(v) Analysis-Test Correlation 

Cantilevered wing models were used in all three tests. It was judged that once the 
wing-winglet interaction was adequately represented, the effect of body and empennage 
on flutter could be accounted for. The pressure model test was designed to collect 


aerodynamic data for both loads and flutter analysis. PiKuie 1 shows the model 
Installation In tlic Boeing Transonic Wind Tunnel (BTWT). P:ess’ir< data was collected for 
a Mach nu.aber-angIo of attack grid for the followrtng configurations; 


A. (1) Clean vvtig width nominal tip 

(11) Clean wing width widnglet at 20* can. an^^le (outboard relative to the 
vertical) 

B. (1) Win(. width nacelle and nominal tip 

(U) Wing width nacelle and (a) Winglet at 20* cant 

(b) Winglet at 10* cant 

(c) Winglet at 0* cant 

C. Configurations described under B above but with the wing sweep angle 
increased by 5* 


Configurations described under B but with the wing sweep angle decreases by 5* 


The pressure data was reduced to sectional data. The wing sectional data was linearized 

width respect to angle of attack to obtain C , and corrected to remove the effect of the 

Tla 

model wing flexibility. The widng sectional data was also linearized with respect tc the 

wwtag sweep angle to obtain C , but was not corrected for the model flexibility effects. 

np 

The widnglet sectional data was similarly linearized without being corrected for the model 
flexibility. The linearized sectional data was used in the flutter arralysis. 

The choice of flutter test configurations and parameters was dictated by the task 
defirdtlon. viz., to develop flutter methodology. Therefore, the test was planned to obtain 
different kinds of flutter modes so that the winglet mass arul aerodynamic effects could 
be separately Identified for each of the flutter modes. The low- speed flutter test was 
designed writh ^ larger number and a wider range of parameters taking advantage of the 
relative ease cf atmospheric low- speed flutter testing compared to high-speed testing. 
The high-speed flutter test was designed after establishing analysis- test correlation for 
the low- speed flutter test. Based on the knowledge derived from the low- speed flutter 
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te«t, a reduced number of configurations arid parameters were selected for testing in the 
high-speed tunnel. The low- speed flutter test was conducted at the General Dynamics, 
Convair Division, San Diego wind tunnel facility. The trzmsonic test was conducted in the 
NASA Langley 16* Transonic Dynamics Tunnel (TDT). A schematic diagram of the wing 
and the wing tips tested, is shown In figure 2. 

The low- speed model wing was of conventional, single-spar construction with wing 
sections perpendicular to the spar. The configurations for the low-speed flutter model 
test were: 

A. (i) Clear, wing (without nacelle) 

(li) Wing with winglet (without nacelle) 

(iii) Winf with winglet mass 

simulator (without nacelle) 

B. (1) Wing with nacelle 

(ii) Wing with nacelle and winglet 

(iii) Wing with nacelle and winglet mass simulator 

C. (i) Wing with nacelle boom 

(li) Wing with nacelle boom and winglet 

(iii) Wing with nacelle boom and winglet mass simulator 

The winglet mass simulator was designed to represent winglet weight, center of gravity 
and inertia properties to help separate winglet inertia and aerodynamic effects. The 
results from configurations with nacelle boom were not used due to good correlation 
obtained for the configurations with nacelle. 

The parameters varied were: 


a. 


c. 


angle of attack, 
model yaw angle, 

wing fuel {0%. 50%, 75%, and 100%), 



f 



d. nacelle strut side bending frequency, 

e. nacelle strut vertical bending frequency, 

f. winglet/simulator cant angle (0®, 10®, 20® relative to the vertical), and 

g. winglet/simulator stiffness. 

The variation of angle of attack and yaw angle was included to evaliiate the static-lift 
effects. The effect of nacelle side bending frequency was found to bo small for the test 
configuration, and is not discussed further in this document. 

The main objective of flutter testing in the NASA Lamgley TDT weis to determine the 
effects of Mach number on flutter characteristics. However, the flutter points obtained 
in a variable density, transonic tunnel depend upon the mass-density ratio as well as the 
Mach effects. Therefore the low-speed model was retested in TDT to determine altitude 
or mass-density ratio effects at lowspeeds. Only two configurations, empty wing with 
nominal nacelle and with and without winglet, were tested. The analysis had shown a 
switch in flutter mode, from rmcelle vertical bending to second wing bending, due to 
decrease in the mass-density ratio. To obtain the mode change in the tunnel, 
mass-density ratio was vauied by testing the configuration wdth winglet in both air and 
freon. The strategy was to show that the mass-density ratio effects, for a winglet 
configured wing, could be predicted at low Mach numbers. The flutter correlation at 
higher Mach numbers could then be evaluated on the basis of compressibility and transonic 
effects. The high-speed model was tested in freon for a Mach range of about 0.6 to 0.91 
and dynamic pressures up to 200 psf . 

The high-speed model was constructed primarily of fibergleiss sandwich components with 
ribs, spars, stringers and skin representing a modem transport vdng. Wing fuel was 
simulated by water. The model was instrumented with 20 accelerometers. 23 pressure 
transducers in two chordwise arrays, and strain gages to monitor wing amd winglet loads. 
The following configurations were selected for testing: 

A. Wing with nacelle and nominal tip 

B. Wing with nacelle and ballasted tip 

C. Wing with nacelle and winglet 
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The ballasted tip configuration W 2 is selected to determine the effect of winglet weight 
separately from winglet aerodynamics. A winglet mass simulator similar to tliat used on 
the low-speed model, would have introdxiced unknown aerodynamic effects at high 
speeds. Therefore, tlie ballast weight v/as incorporated inside the vdng contour resulting 
in a wing tip aerodynamically identical to the nominal tip. The test par 2 uneters selected 
were: 

a. wing fuel (empty and f\ill), 

b. nacelle strut vertical bending frequency, 

c. winglet cant angle (0® and 20* relative to the vertical), and 

d. angle of attack. 

Two nacelle strut vertical bending springs were used. The nominal strut vertical bending 
spring (nominal nacelle) and the softer strut vertical bending spring (soft nacelle) gave 
rise to different flutter characteristics due to differences in coupling of nacelle motion 
with inboard wing torsion. A series of high angle of attack runs, within the model load 
limits, was run to verify that there were no single-degree-of-freedom instabilities at 
transonic speeds. 

This volume pertains to the flutter test conducted in the NASA Langley tunnel. The 
highlights are covered in the main body. Appendices A, B and C contain sufficient data, 
in the form of figures and tables, to allow an independent analysis. Appendix D contains 
procedure used to modify calculated stiffness matrix. A summary of experimental results 
is tabulated in Appendix E. 


3.0 Mass- Density Ratio Effects at Low Mach Numbers 

The low-speed model tested earlier in the Convair tunnel, was retested in the TDT. Only 
two configurations, empty wing with nominal nacelle and with and without winglet, were 
tested. The configurations tested were identical to the similar configurations tested 
earlier at Convair e.:cept for the following; 
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(a) The model was wall mounted and supported from the balance in the TDT. 

The body fairing of the high-speed model was used. In the Convair tirmel. the 
model was mounted on a stiffened body supported by an A-frame bolted to the 
tunnel floor at the centerline of the tunnel. The model test frequencies did not 
change significantly between the two installations. Table 1 lists the two sets of 
frequencies and the analyatical frequencies for the configuration with winglet. 
Figure 3 shows a photograph of the low- speed model installed in TDT. 

(b) The model tested in the Convair tunnel had shims installed in the wmg sections 
to simulate the wing twist distribution for most of the runs. The model was 
installed in the TDT. without any shims. The effects of wing twist for the 
winglet configuration with 75% fuel ceise was evaluated in the Convair test and 
is summarized below; 


FLUTTER 


TUNNEL 

WING 

TIP 

SHIMS 

SPEED 

(KTAS) 

DYN PR 
(PSF) 

FREQ 

(Hz) 

CONVAIR 

WINGLET 

YES 

90.9 

27.4 

8.7 

CONVAIR 

WINGLET 

NO 

88.8 

25.8 

8.6 

CONVAIR 

NOMINAL 

YES 

97.9 

31.7 

8.7 

CONVAIR 

NOMINAL 

NO 

96.6 

30.8 

8.6 


Since the repeatability of flutter speed was determined tc be within 1 KTAS, there 
appears to be a small drop in the flutter speed due to removal of the shims from the wing 
sections for the configuration with winglet. The repeatability of flutter speeds between 
the Convair and TDT (air) for the two empty fuel configurations is summarized below: 





DENSITY 





WING 


(SLUGS/ 

SPEED 

FLUTTER FREQ 

TUNNEL 

TIP 

SHIMS 

FT3)x 1o3 

(KTAS) 

DYN PRESS(Hz) 
(Ib/ft*^) 


CONVAIR 

WINGLET 

YES 

2.329 

89.2 

26.4 

8.8 

TDT 

WINGLET 

NO 

2.309 

86.0 

24.3 

9.0 

CONVAIR 

NOMINAL 

YES 

2.322 

96.5 

30.8 

8.7 

TDT 

NOMINAL 

NO 

2.349 

91.2 

27.8 

9.0 


The difference in the flutter speed for the configuration without winglet is about 5%, and 
is higher than the configuration with winglet. The higher difference in flutter speeds for 
configuration with nominal tip was not looked into in detail as more emphasis was placed 
on the configuration with winglet. The mass- density ratio effects were obtained and 
analyzed in more detail for the configuration with winglet as described below. 

The analysis had shown a switch in flutter mode, from nacelle vertical bending to second 
wing bending, due to decrease in the mass-density ratio. The range of mass-density ratio, 
to affect the mode change in the tunnel, was achieved by testing the configuration with 
winglet in both air and freon. Figure 4a shows the analysis- test correlation as a function 
of mass-density ratio for all test points. In Figure 4b, the data points are shown for 
mass-density ratios up to 50 to show more cleawly the switching of flutter modes. The 
analysis is able to predict the trend correctly auid shows good correlation with the test 
results. The switch in the flutter mode occurred at higher mass-density ratio in the test 
than shown by analysis. A small difference in actual and predicted damping could explain 
this difference. The amalytical results were calculated using the post- test model 
described in Section 6 of Volume I except that the analysis frequencies were adjusted to 
match the model GVT. This auljustmerit is equivalent to about 1.5 KTAS increase in 
anal 3 dical flutter speed. It was concluded that the mass-density ratio effects can be 
predicted with acceptable accuracy for winglet configured wings at low Mach numbers. 



/ 
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4.0 Description of High-Speed Test 


The model was installed in the NASA Langley TDT as shown in Figure 5. The model was 
supported on the NASA balance mounted on the wall turntable. A badlasted wing tip was 
fabricated such that it could replace the nominal wing tip. The weight and Body Station 
coordinate of center-of-gravity location of the ballasted tip were similar to the winglet: 


WEIGHT CG LOCATION (INCHES IN WRP 

(LBS) RELATIVE TO LEADING EDGE OF WING TIP) 

AFT OUTBOARD 


BALLASTED TIP 

0.350 

5.00 

0.40 

WINGLET 

0.378 

5.36 

1.17 

NOMINAL TIP 

0.0198 

3.96 

0.30 


The nominal and badlasted tips were aerodynamically identical to each other. Each one of 
the three tips could be attached to two hard points m the wing tip structure . 

The nominal nacelle configuration was a strut - spring combination corresponding to 
cantilevered nacelle vertical bending frequency of 24.7 Hz. The soft nacelle configuration 
was a strut - spring combination corresponding to a cantilevered nacelle vertical bending 
frequency of 15.99 Hz. The nominal winglet configuration was a 20° cant winglet with a 
cantilevered winglet frequency of 93.0 Hz. 

The test procedure adopted reflects the emphasis on flut«.er correlation rather than 
flutter clearance. A comprehensive model GVT was conducted in the Boeing Structures 
Dynamics Laboratory (SDL) prior to the wind tunnel test. The model frequencies were 
also determined for the tunnel installation. The model responses were monitored during 
the test. Some of the highlights of the test procedure are described below. 
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The cantilevered model configurations tested in the SDL for mode shapes were: 




BODY 

WING 




bdPPORT 

SHELL 

FUEL 

NACELLE 

WINGLET 

a) 

STEEL PLATE 

OFF 

EMPTY 

OFF 

OFF 

b) 

NASA BALANCE 

OFF 

EMPTY 

OFF 

OFF 

c) 

NASA BALANCE 

ON 

EMPTY 

OFF 

OFF 

d) 

NASA BALANCE 

ON 

EMPTY 

NOMINAL 

OFF 

e) 

NASA BALANCE 

ON 

EMPTY 

NOMINAL 

NOMINAL 

0 

NASA BALANCE 

ON 

FULL 

NOMINAL 

OFF 

g) 

NASA BALANCE 

ON 

FULL 

NOMINAL 

NOMINAL 

h) 

NASA BALANCE 

ON 

EMPTY 

OFF 

NOMINAL 


The steel plate as well as the NASA Balance were supported from a strongback. A 
hammer test for wing-nacelle-winglet configuration was conducted in the SDL and it was 
decided that the hammer test, instead of shaker, Aill oe used in the tunnel to obtain 
model frequencies in still air. An instrumented hammer wais used for exciting the model 
to obtain its frequency response. The frequency spectnim of each configuration in the 
tunnel was examined to verify that the model was not damaged amd the tunnel installation 
was proper. This procedure allowed detection of anomalies in the model and helped in 
isolation and correction of the cause of any differences. 

Figure 6 shows the tunnel characteristics for freon operation. The total pressure (H) 
curves roughly correspond to constant stagnation density lines in the tunnel. The most 
efficient tunnel operation is achieved in a tunnel "run" by operating the tunnel alor.g 
constant H curves which results in a simultaneous incre£u>e in the Mach numiber and 
dynamic pressure. A run was terminated if one of the following four conditions was 
reached: (i) Mach 0.91, or (li) dynamic pressure of 200 psf, or (iii) excessive model 
response amplitudes either due to buffet or low damping, or (iv) onset of flutter. The 
maximum Mach number, dynamic pressure and other tunnel parameters well as the 
frequencies of significant responses and reason for terminating the run were recorded. 
This procedure was followed throughout the test. The results of pre-test analysis were 
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used to select the tunnel runs. For a typical predicted flutter boundary example, also 
^own in Figure 6, runs might be made along H - 300, 400 and 600 psf. If the test flutter 
points obtained were judged to correlate with the prediction, no further runs were 
regarded as necessary for that particular configuration; otherwise additional nins were 
made. A decision was made to not spend the tunnel time in precisely defining the 
transonic bucket. Instead, the limited tunnel time was utilized to test as mmy 
configurations as possible. This strategy proved to be successful based on the number of 
configurations tested and flutter points obtained in the tunnel. 

During each nm the model responses were monitored by strip chart traces, and reduced to 
power spectra, "damping Indicator", and cascade plots. The acceleration power spectrum 
of either the wing tip vertical and fore- aft accelerometers (or both) were displayed in 
real time with updates every second. The inverse of the amplitude of the highest peak of 
the power spectrum, for wing tip vertical accelerometer, was plotted, also in real time, 
versus Mach number. This was called a "damping indicator" or "relative damping" plot. 
The cascade power spectra were plotted for almost every run and were available within a 
few minutes after the run. 

Two DRAS (Dynamic Response Actuated Switch) units were employed to safeguard 
against excessive model acceleration amplitudes. One of the DRAS units was set to 
actuate opening of the four tunnel by-pass valves for quick shutdown at a preset, 
sustained amplitude. The second unit was set up to switch a red warning light at a fixed 
percent of the shut-off amplitude. The hook up of the DRAS units was initially permitted 
by NASA on a trial basis. There was a concern about DRAS unit repeatedly shutting-off 
the tunnel prematurely. The system worked well in practice and the DRAS unit remained 
in the shut-off loop throughout the test. 


5.0 Analytical Representation 

The analytical representation used W 2 is Identical to the low- speed model (Volume 1). The 
built- up, high-speed model wing was structurally represented by finite beam elements 
(elastic axis) as if the wing were of single- spar cc struction. The nacelle and strut were 


attached as rigid, lump-masses to the wing elastic axis. The winglet and ballasted tip 
were represented as separate substructures using branch mode representation. The 
cantilevered nacelle strut and winglet test frequencies and mode shapes were input as 
assumed modes. The mass distributions of the model wing and winglet were calculated, 
and the total mass and inertias were individually verified with the measured values. The 
nominal and ballasted tip mass and inertia properties were measured. The calibrated 
model stiffness properties were used to improve correlation with the results of the model 
GVT. This data is included in Appendix A to allow independent analysis. 

The aerodynamic representation for flutter analysis was based on the strip-theory 
aerodynamics (AFl program-Ref 8). The sectional aerodynamics data was derived, from 
two wind tuimel tests, for Mach numbers of 0.4, 0.65, 0.80, 0.88, and 0.91. There were 
minor differences between the pressure model and the high-speed flutter model. 
Therefore, an earlier wind tunnel test for wing-nacelle cordTguration was used as a basis 
for sectional aerodynamics data. To obtain the sectional data for configurations with 
winglet, the difference due to winglet from the later test was algebraically added to the 
sectional data from the earlier test. The sectional data for the five selected Mach 
numbers is included in Appendix B. In order to get a theoretical sectional aerodynamic 

data. DUBLAT (doublet lattice program - from Ref 8) was used for steady flow. The 

c 

theoretical sectional data obtained are also included in Appendix B. The nacelle n ^ and 
r 

y p values used were 0.052 2 ind 0.042 respectively at M = .4 and changed very little at 
high Mach number. 

The flutter solutions were obtained at Mach numbers of 0.4, 0.65, 0.80, 0.88, and 0.91 and 
five densities for each configuration. The Qutter dynamic pressures were plotted on 
tunnel charts. The match-point solutions were determined, for each flutter mode, 
corresponding to structural damping (g) of 0.0 and 0.03. 


6.0 Correlation With Model GVT Results 

The model GVT results were used to modify the analytical model to improve correlation 
with the test mode shapes and frequencies. The GVT results for the clean wing (without 
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nacelle and with nominal tip) configuration were used to modify the analytical 
representation. Assuming the analytical mode shapes to exactly match the test mode 
shapes for the clean wing, the analytical stiffness matrix was modified based on matching 
the frequencies. Table 2 lists the frequencies for the clean wing. The node line plots for 
mode shapes from the modified analysis and the test are included in Appendix C for clean 
wing, wing-nacelle, and wing-nacelle-winglet configurations. The procedure for 
modifying the stiffness matrix is described in Appendix D. 

The modified stiffness matrix was used for all configurations. This was possible because 
of the modeling approach described in Section 5. Tables 3 and 4 list the frequencies for 
wing-nacelle, wing-nacelle-ballasted tip and wing-nacelle-winglet configurations for 
empty and full wing fuel. The differences in GVT frequencies between model installation 
in the Structures Dynamics Laboratory (SDL) and the tunnel are attributed to the tunnel 
turntable. The frequencies for empty wing configurations with soft nacelle are shown in 
Table 5. 

The modal correlation between the analysis and GVT was considered to be reasonably good 
specially when considering the difficulties involved in stiffness calibration. The only 
significant difference in the modal correlation was found to be for the wing (full)- nacelle 
(nominal)- winglet (20 deg) configuration for the "chordwise" mode. For this mode, the 
wing fore-aft motion is dominant with significant coupling with outboard wing bending 
and torsion. The analytical frequency is about 2 Hz nigher than the frequency of 20.7 Hz 
obtained during tunnel GVT. In the post- test analysis, this difference was found to be 
significant for some flutter modes. This is further discussed in Section 9. 


7.0 Flutter Test Results and Correlation 

The test results obtained are shown in figures 7a to 7c. There was a significant reduction 
in flutter dynamic pressure (Qp) due to the winglet aerodynamic effects. For the 
configuration with nominal nacelle strut and empty fuel (fig. 7b). the effect of the 
ballasted tip was to slightly lower the flutter boundary except at higher Mach numbers. 
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However, for the configuration with the nominal nacelle strut and full fuel (fig. 7c). the 
effect of the ballasted tip was to cause a low-damped mode to occur at a slightly lower 
dynamic pressure. The reduction in Qp due to winglet aerodynamic effects was more 
pronoxmced for this case. The configuration with the soft nacelle strut and empty fuel 
(fig. 7d) showed trends similar to the configuration with nominal nacelle strut. The effect 
of winglet cant angle shown in figure 7e. was foimd to be similar to that for the low-speed 
model. The differences in the effects of winglet aerodynamics on different configurations 
were primarily due to the flutter modes. The four flutter modes encountered were similar 
to the four flutter mechanisms foimd for the low-speed model. The flutter test results 
are summarized in a tabular form in Appendix E. An angle-o^ -attack variation series was 
run, within model load limits, over a Q-M range representative of scaled flight envelope 
for the wing (empty)- nacelle (nominal)-wlnglet (20 deg) configuration. The model load 
limits were -80 lbs to 180 lbs. No single degree-of-freedom instability was found to exist. 

The analysis-test correlation obtained is shown in figures 8-11. A short discussion related 
to each configuration is presented below. 

(a) Wing (Empty)- Nacelle (Nominal) configuration results are shown in Figure 8a. 
The nacelle vertical bending mode was found to flutter. The analysis also 
predicts the flutter to occur in the nacelle vertical bending mode. The 
analytical Op-M flutter boundary appears to have similar shape as the test, but 
the analysis is conservative. 

(b) Wing (Empty)-Nacelle (Nominal)- Ballasted Tip configuration results are shown 
In Figure 8b. The flutter stUl occurs in the nacelle vertical bending mode. 
However, the analysis is slightly unconservative rather than conservative as for 
the nominal tip configuration. 

(c) Wing (Empty)-Nacelle (Nominal) -Winglet (Nominal) configuration results are 
shown in Figure 8c. At the two higher Mach numbers, M » 0.77 and .828, the 
flutter occurred in nacelle vertical bending mode. At Mach .66, the model 
response showed high amplitude in 17.6 Hz nacelle vertical bending mode and 
22.3 Hz second wing bending mode. 
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The ratio of acceleration amplitude squared, of the 22.3 Hz to the 17.6 Hz 
mode is 1.37. This ratio is based on a spectrum derived from exponential 
averaging, with overlap processing, of ten ensembles of five seconds each. The 
corresponding ratio of displacement amplitudes is 0.72. Therefore, it is possible 
to classify the flutter mode as a second wing bending mode based on 
acceleration response, or as a nacelle vertical bending mode based on 
displacement response. However, since acceleration response is generally used 
in flutter testing, the flutter mode was designated as second wing bending 
mode. The analysis-test correlation for nacelle vertical bending mode is 
satisfactory. The correlation at M = 0.66 is also satisfactory, since the test 
point is interpreted as being a combination of nacelle vertical bending second 
wing bending response. 

(d) Wing (Full)- Nacelle (Nominal) configuration results are shown in Figure 9a. The 
wing chordwise bending mode disappears at higher damping (g). Two runs were 
made and one flutter point for wing tip mode was obtained. The analysis- test 
correlation is satisfactory. 

(e) Wing (Full)- Nacelle (Nominal)- Ballasted Tip configuration results are shown in 
Figure 9b. Two runs were made. No flutter points were obtained although the 
pass at higlier dynamic pressure resulted in some low damped response in second 
wing bending mode. 

(f) Wing (Full)- Nacelle (Nominal)- Winglet (Nominal) configuration results are 
shown In Figure 9c. This is the most complicated configuration in terms of 
sorting out the flutter modes. Three flutter modes (nacelle vertical bending, 
wing chordwise bending and wing tip) were observed. At Mach .856, there was 
distinct beating between the 18.5 Hz (nacelle vertical bending) and the 19.1 Hz 
(wing chordwise bending) modes. At Macn .79, response in both these modes is 
apparent. The higher frequency wing tip mode was observed for the test points 
at ‘lach .73 and .644. The analytical results match fairly well for the nacelle 
vertical bending and v/ing chordwise mode at the two higher Mach numbers. 
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However, the analysis appears too conservative for the wing chordwlse mode 
aivd unconservative for the wing tip mode. Considerable aiulytical effort was 
devoted In understanding the sensitivities of the wing chordwlse and tip modes. 
The results are discxissed in Section 9. 

(g) Wing (Empty)-Nacelle (Soft) configuration results are shown m Figure 10a. At 
the higher Mach numbers, low-damped response at several frequencies was 
observed. The flutter obtained was in second wing bending mode. The 
analysis- test correlation is satisfactory. The area included in the g > 0 
boundary for the nacelle vertical bending mode is the region of Instability for 
that mode. 

(h) Wing (Empty)-NAcelle (Soft)-Ballasted Tip configuration results are shown in 
Figure 10b. Tlie flutter characteristics and the nature of correlation is similar 
to case (g) above. 

(1) Wing (Empty)-Naoelle (Soft)-Winglet (Nominal) configuration results are shown 
in Figure 10c. The flutter speeds are lower compared to cases (g) emd (h) 
above. The analysis-test correlation is satisfactory. An interesting feature is 
that both anal 3 Tsi? and test, show presence of the nacelle vertical bending and 
second wing bending modes in close proximity to each other. 

(j) Wing (Empty)- Nacelle (Nominal)- Winglet (O’ cant) configuration results are 
shown in Figure 11. The flutter characteristics did not show a significant 
difference due to change of winglet cant angle from 20* to 0*. 

The effect of the winglet was to reduce the flutter dynamic pressure. In order to get a 
quantitative effect, the analytical flutter results for wing (empty)- nacelle (nominal) with 
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the three wing tips are shown in Figures 12a to 12c for a density of 1.11 X 10 slugs/ cu 

ft. For the nacelle vertical bending mode at M » .88, the effect of the tip weight is to 
increase the flutter dynamic pressure (g - .03) by 7% relative to the nominal tip 
configuration. The aerodynamic effect of winglet is to reduce the flutter dynamic 
pressure ( g » .03) by 14% relative to the ballasted tip configuration. Thus the net 


reduction In the analytical flutter dynamic pressure due to the combined effect of winglet 
weight and aerodynamics Is relatively small for the empty fuel configuration. For 
different wing configurations, different fuel conditions or different modes, the effect of 
the winglet welgjjt and aerodynamics will obviously be different. 


8.0 Reduction of Test Data 

The data obtained during the test was critically evaluated after the test. The reason for 
stopping each run was reviewed by examiiUng the strip charts and cascade plots. Also 
post- test data reduction was done for selected runs. The post- test data reduction 
consisted of: 
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(a) Plotting of calibrated l/A (from power spectrums) versus Mach number for 
maximum amplitude in each of the three. r*lected frequency bands. 

(b) Calibrated time histories for eight seconds near end of tn«. nm with the 
maximum response near the middle of the eight seconds. 

(c) Flutter modes were derived from Fourier analysis of one second of maximum 
response described in (b). The wing tip response was used as a reference in 
defining the phase relatior\s. 

Figures 13, 14, ana 15 present examoles of the information provided by (a), (b), and (c), 
respectively. The following observations were made from review of the test data. 

(a) The flutter dynamic pressure and Mach number recorded In the tunnel did not 
warrant much of an adjustment. The "damping indicator" versus Mach number 

plots could be used to extrapolate to a Mach number corresponding to a 

2 2 
selected level of 1/A . It is not feasible to extrapolate to a 1/A » 0 

corresponding to flutter because of (i) some scatter in the "damping Indicator" 

versus Mach plots, (li) the accuracy of .01 associated with determination of the 

tunnel Mach number and (Hi) the accuracy associated with reading the recorded 



analog voltage signal corresponding to the Mach number. It was confirmed that the 
tunnel was stopped at about the same level of acceleration response for most of the 
flutter runs. Therefore, no modification was made to flutter Mach numbers and 
dynamic pressures recorded during the test. 

(b) Comparisons of time traces a^^p^oaching flutter for runs where the response 
frequencies were significantly different, do not indicate significant differences 
in phase relationsMp c^ one accelerometer to the other. From the low-speed 
model responses, it was possible to identify different phase characteristics for 
different flutter modes. This is much more difficult to d.'' from the high-speed 
model responses. Possible reasons could be that (i) the higher turbulence 
associated with higher d 5 mamic pressures excites many of the lower frequency 
modes, and (ii) the sustained oscillations are not maintained for a sufficient 
time to clejun up the response. Therefore, the' response frequency remains the 
primary means of identifying the flutter mode. 


9.0 Analytical Seiisitivity Studies 

An analytical sensitivity study was conducted to evaluate the effect of selected 
parameters on analysis- test correlation. The primary configuration for the sensitivity 
study was the wing (full)-nacelle (nominal)- winglet (nominal) configuration. 'This was 
judged to be the most interesting configuration tested since three flutter mecli«irilsnis 
were observed. In addition, this configuration was found to be sensitive to the 
characteristics of the wing chordwise bending mode as described later in this section. 

The following parametric variations resulted in small changes to the flutter results, and 
were judged to be not significant: 

(a) Win^ elastic axis location was varied as shown in Figure 16. 
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(b) Number of aerodynamic strips in the AFl program were varied, 

(c) The static-lift effect was included in the flutter analysis using SLOAEF 
program. 

(d) The modification of stiffness matrix (Appendix D). 

The flutter speeds were found to be sensitive to the following parameters: 

(a) Structural - Wing chordwise bending, and 

(b) aerodynamic - spanwise distribution of static lift-curve slope 

and aerodynamic center. 

The sensitivity to wing chordwise bending was considered to be somewhat unusual, and is 
believed to be brought about by the combination of winglet and fuel. The importance of 
the wing chordvdse bending mode can be seen from plots of flutter dynamic pressure 
versus wing chordwise bending frequency for Mach (density) combinations of 0.65 (3.50 X 
10”^ slugs/cu ft.) and 0.88 (1.11 X 10'^ slugs/cu ft.) in Figures 17a and 17b. The wing 
chordwise beruilng mode has a significant wing tip vertical motion component which 
accounts for its effect on flutter speeds. 

The chordwise bending stiffness was modified (see Figure A2) to evaluate the effect of 
change in stiffness distribution. The modified stiffness probably was a better 
representation of the model. The chordwise bending stiffness hats significant effect on the 
frequencies of wing chordwise and torsion modes. The resulting wing chordwise bending 
and torsion frequencies for the wing (full)- nacelle (nominal)- winglet (nominal) were 20.3 
Hz and 40.3 Hz. respectively. The corresponding frequencies for the reference analysis 
were 22.87 Hz and 41.89 Hz (Figure C5). The corresponding test frequencies were 20.7 Hz 
and 42.8 Hz (TDT Installation), respectively. Thus the modification to the wing chordwise 
bending stiffness Improved the GVT frequency correlation for the chordwise mode with 
some deterioration for the torsion mode. The sensitivity of flutter dynamic pressure to 
chordwise bending stiffness, along with other parameters, is discussed below. 
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The flutter dynamic presjwe (Qp) at selected Mach number - density combinations is 
compared in Taole 6. Four Mach niimber - density combinations were selected mainly 
based on their proximity to the test points obtained. These were Mach = 0.4, 0.65, 0.80, 
and 0.88 and corresponding densities of 1.11 X 10 3.50 X 10 1.50 X 10~^ and 1.11 X 

_3 

10 slugs/cu ft, respectively. The Qp for reference analysis is tabulated for comparison. 
The base for sensitivity analysis is different from the reference analysis in the chordwise 
bending stiffness. 

Variation 1 shows the effect of "tuned” frequencies. The analytical frequencies were 

"tuned" to match the GVT frequencies. The eflect is primarily due to the wing torsion 

mode frequency change from 40.3 Hz to 42.8 Hz, and results in increasing Qp of the tip 

mode. In Variation 2, there is a 10% increase in wing C for i) - .538 to 1.0. There is a 

drop in Qp for all three modes. The effect of shift in wing aerodynamic center (tj= .538 to 

1.0) by .05C arid .IOC forward, is also seen to be significant from results tabulated under 

Variations 3 and 4. As expected, Q„ drops for all three modes. The stiffness Variation 5, 

r 

shows the effect of including the stiffness matrix modification based on the cantilevered 
wing only frequencies. There appears to be a further drop in Qp for the nacelle vertical 
bending and wing chordwise bending modes. This requires further evaluation to understand 
the reasons for the significant effect. 

The effect of using doublet-lattice aerodynanucs program (DUBLAT) rather than the strip 
theory program (AFl) was evaluated. No empirical corrections were used. The results 
tabulated under Variation 6 should be compared to Variation 2 for M = .4 and .65, and to 
Variation I for M » .8. The reason being that the test C sectional distributions used for 

da 

M = .4 and .65 are about 10% lower than the corresponding theoretical DUBLAT 
distributions. For M .» .8, the test and theoretical C distributions are similar. It is not 
surprising that the flutter dynamic pressure predicted by DUBLAT at M = .88 is 
significantly higher than the AFl results as well as the test results. The DUBLAT resxilts 
at Mach 0.4 and 0.65 appeau* to be in the right range. However, the reason for DUBLAT 
predicted flutter dyii .inic pressure at M ^ .80 being significantly higher than the results 
from the base run, needs to be investigated. 
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The effect of chordwise mode shape appears to be significant based on preliminary 
assessment. As described earlier, there is a significant wing bending and torsion motion 
associated with the wing chordwise bending mode for the wing (full)-nacelle 
(nominal)-winglet (nominal) configuration. An attempt to use experimentally measured 
mode shape was initiated, but has not been completed. It was found that the wing twist 
could not be reliably reduced from the measured data as it was sensitive to small changes 
in accelerometer readings. The error bounds for the accelerometers are not known with 
sufficient accuracy to enable evaluation of the quality of wing twist information obtained 
from the displacement data. 

The application of compressibility correction or ”C^ correction" was compared to direct 
solution using API with empirical sectional data for appropriate Mach numbers. The 
flutter at any Mach number. ( n\ay be determined as 


Q 





] 


2 

where C is generally determined from the wind tunnel test data, and M, is selected to be 

Incompressible Mach number. It has been customary to use M^ = 0.4. Figures 18a to 18d 

show comparisons of two methods of solutions for wing-nacelle and wing- nacelle- winglet 

configurations at a density of 1.11 X 10 slugs/cu ft. The comparison is shown for two 

flutter modes, nacelle vertical bending and second wing bending mode. For the 

wing nacelle configuration, the C correction resulted in higher Q— at transonic Mach 

c r 

numbers for both modes. However, the nacelle mode is softer at higher Mach numbers 

using actual Mach solution compared to the solution obtained with C correction. For the 

c 

wing- nacelle-winglet configuration, the C solution approximated fairly well the actual 

c 

Mach solution. These comparisons have been made for specific configurations and 

altitude. No general conclusiot^s aure warranted except one. The "C correction" may give 

c 
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results different from the actual Mach solution, and may not be always conservative. 
However, the simple approach may be useful for preliminary evaluation of test 
configurations for the purpose of planning the test. 

10. Conclusions and Recommendations 

The test program has been successful in creating the data bsise for flutter characteristics 
of vdnglet configured wing for a twin-engine configuration. The four flutter mechanisms 
predicted by analysis were obtained in the tunnel. The number of flutter test points for the 
ten high-speed configurations and two low-speed configurations, obtained in the tunnel 
cover a wide range of altitudes and Mach numbers. This provides an excellent reference 
for evaluation of analytical correlation for a configuration with and vdthout winglet. 

The mass-density ratio effects at low Mach numbers were correlated (auialysis vs. test) 
satisfactorily over a wide range. The application of conventional analysis proved to be 
satisfactory through the transonic Mach regime. It was not surprising to find that 
theoretical doublet-lattice analysis gave unconservative answers at M - .8 and .88. It was 
concluded that the flutter characteristics of a winglet configured high aspect ratio wing 
can be satisfactorily predicted with careful application of existing methods for a 
twin-engine airplane configuration. The wing chordwise bending mode for certain 
configurations can be expected to have significamt wing bending and torsion motion. It is 
indicated by the present study that this coupling effect is important. 

It is recommended that the experirriental and analytical data base established in this 
program be used to awivantage. The number and diversity of flutter test points and the 
correlation established with simple methods, should be used to evaluate state-of-the-art 
transonic codes. There are many examples published where two or three 
degrees- of-freedom systems or simplified representations have been studied for transonic 
effects using very expensive codes. It is believed that the time has come to make a real 
effort using the data from a realistic configuration to determine the advantages and costs 
of applying transonic codes. 
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TABLE 1. LOW-SPEED MODEL FREQUENCIES 
EFFECT DUE TO TUNNEL INSTALLATION 
WING(EMPTY) - NACELLE (NOMINAL) - WINGLET (20 DEG) 

MODEL FREQUENCIES (Hz) 
FOR INSTALLATION IN 


MODE 

CONVAIR TUNNEL 

NASA TDT 

ANALYSIS 

1st WING BENDING 

3.93 

4.0 

3.78 

NACELLE SIDE BENDING 

- 

8.0 

7.83 

NACELLE VERTICAL BENDING 

9.45 

9.5 

9.13 

WING CHORDWISE BENDING 

11.74 

13.49 

12.4 

2nd WING BENDING 

12.19 

13.0 

11.93 

NACELLE ROLL 

17.88 

- 

- 

1st WING CHORDWISE BENDING + 
0/B WING TORSION 

23.46 

24.6 

22.85 

0/B WING TORSION + WING 
CHORDWISE BENDING 

23.80 

24.17 

25.6 

1st WING TORSION 

27.00 

26.6 

31.73 
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TABLE 2: HIGH SPEED MODEL. CORRELATION OF ANALYSIS AND TEST VIBRATION 
FREQUENCIES (HZ) FOR CLEAN WING (EMPTY) 


MODE TEST ANALYSIS 



RIGID 

PLATE 

1st Wing Bending 

7.81 

2nd Wing Bending 

25.00 

1st Wing Chordwise Bending 

34.00 

3rd Wing Bending 

52.68 

1st Wing Torsion 

58.42 

Higher Mode 

86.7 

Higher Mode 

96.10 


NASA NASA WITH 

BALANCE BALANCE STIFFNESS MOD 


7.80 

7.61 

7.80 

24.70 

23.93 

24.70 

32 02 

34.90 

32.02 

52.12 

52.61 

52.12 

58.08 

57.70 

58.08 

85.88 

90.83 

85.88 

94.12 

95.91 

94.12 
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TABLE 6 - SUMMARY OF SENSITIVITY STUDIES. WING(FULI)-NACELLE(N0M1NAL)-WINGLET(20°) 
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NACELLE 


aspect ratio 

7*88 

QUARTER CHORD SWEEP 

-31° 

LOW. SPEED MODEL SPAN 
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HIGH-SPEED MODEL SPAN 

93.5” 
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(WINGLET/ EMPTY WING)HEIGHT .012 
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FIG. 2 MODEL WING AND WING TIPS 
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FIGURE 7b EFFECT OF WINGTIP COMFIGUPATIOH OH TEST FLUTTER BOUNDARY, 
WING (EMPTY)- NACELLE (NOMINAL) 
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FIGURE 7c effect OF WINGTIP COHFIGUPATION ON TEST FLUTTER BOUNDARY, 
WING (FULL) - NACELLE (NOMINAL) 
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FIGURE 7d EFFECT OF WIMGTIP CCNFIGUPATIOM OH TEST FLUTTER BOUNDARY, 

WING (EMPTY) - NACELLE (SOFT) 
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FIGURE 7e EFFECT OF WIHGI.ET CANT ANGLE ON TEST FLUTTER BOUNDARY, 
WING (EMPTY) - NACELLE (NOMINAL) - WINGLET 
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FIGURE 3a FLUTTER CORRELATICU FOR WING (EMPTY) - NACELLE (NOMINAL) 
NOMINAL TIR 
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FIGURE 10b FLUTTER CORRELATION FOR WING (EMPTY) - NACELLE (SOFT) - 
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FIGURE 18a COMPRESSIBILITY CORRECTION AND FLUTTER DYNAMIC PRESSURE. 
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ALL INERTIAS ARE ABOUT C.G. 


FIGURE A5 FULL FUEL MASS & INERTIA PROPERTIES 
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CANTILEVER NACELLE FREQUENCIES AND MODE SHAPES * 

No . MODE 

FREQ (HZ) 

MODE SHAPE 

1 NAC SIDE BNDG 

2 NAC VERT BNDG 

3 NAC ROLL 

15.74 

24.70(N0M),15.99(S0FT) 

28.97 

TY 1.0 RX .152 RZ -.111 
TX 1.00 TZ -.968 RY -.179 
TX .72336 TY 1.0 TZ -.89376 
RX 6.591 RY -.181 RZ 2.509 



CANTILEVER WINGLET MODE SHAPES** 
(93.00 HZ) 

WINGLET NODE 
AT WBL 

TZ 

RX 

RY 

1.351 

.003174 

.00061688 

-.0003185 

4.197 

.055950 

.0029368 

-.0015398 

6.760 

.20376 

. 006447 

-.0034839 

9.051 

.45002 

.01098 

-.0063343 

11.939 

1.0000 

.018689 

-.011197 


*IN GLOBAL FRAME 
**IN WINGLET REF FRAME 


FIGURE A6 


CANTILEVERED NACELLE AND WINGLET FREQUENCIES AND MODE SHAPES 



APPENDIX B 


AERODYNAMIC DATA 


WING ~ PER DEG 












BASELINE(HING WITH NACELLE) 
WITH WINGLET(20° CANT) 

WITH W1NGLET(10° CANT) 
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fcTA , PERCENT WING SEMI -SPAN 















i ETA, PERCENT WINGLET SEMI -SPAN 


LEADING 

EDGE 

















FIGURE B15e THEORETICAL WING C DISTRIBUTION, WING-NACELLE-WINGLET(20 DEG.) 

” a 

M - 0.91 
102 

















LEADING EDGE WING + NACELLE + WINGLET(NOM) OUBLAT ANALYSIS 







APPENDIX C 


VIBRATION FREQUENCIES & NODE LINES 



i 


i 1 


WING BENDING 



TEST (HZ.) 7.81 
ANALYSIS (HZ.) 7.80 


WING BENDING 



24.70 


1^^ WING CHORDWISE 



32.02 


WING BENDING 



TEST (HZ.) 52.68 
ANALYSIS (HZ.) 52.12 


HIGHER MODE 



58.42 

58.08 


HIGHER MODE 



85.88 


HIGHER MODE 



TEST(HZ.) 96.10 
ANALYSIS (KZ.) 94.12 


LEGEND; 

NODE LINES 

test 

ANALYSIS 

*TEST NODE LINE NOT 
AVAILABLE 


FIGURE Cl MEASURED S CALCULATED FPEQUEUCIES i NODE LINES FOR WING (EMPTY) 
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1st WING BENDING 



TEST(HZ.) 7.72 

ANALYSIS(HZ.) 7.76 


NACELLE SIDE BENDING 



15.26 


2nd WING BENDING 



ANALYSIS(HZ.) 24.15 


NACELLE ROLL 



29.49 

29.80 


3rd WING BENDING 



ANALYSIS(HZ.) 42.76 


1st WING TORSION 



57.19 


NACELLE VERTICAL 



19.82 

18.60 


1st WING CHOROWISE 



30.47 

32.07 


LEGEND: 

NODE LINES 

TEST 

ANALYSIS 


FIGURE r? f'EASUPED I CALCULATED FPECUEHCIES A NODE LINES FOR 
WING (EMPTY) - NACELLE(NPMINAL) 



TESKHZ. 

ANALYSIS(HZ. 


TEST(HZ 

ANALYSIS(HZ 


TEST(HZ. 

ANALYSIS(HZ. 


1st WIK'G BENDING 



) 6.72 

2nd WING BENDING 




) 21.08 
3rd WING BENDING 

A 

/ f 


/ / 



NACELLE SIDE BENDING 



15.23 

15.19 


1st WING CHORDWISE 



26.64 


1st WING TORSION 



47.01 


NACELLE VERTICAL 



18.54 


NACELLE ROLL 



29.98 


LEGEND; 

NODE LINES 

TEST 

ANALYSIS 


FIGURE C3 MEASURED A CALCULATED FREOUENCIES A ,'''CDE LINES FOP 
WING (EMPTY) - NACELLE (NOMINAL)- WINGLET (20 DEG) 
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A 


1st WING BENDING 



TEST(HZ.) 6.06 

ANALYSIS(H2,) 6.06 


NACELLE SIDE BENDING 



14.98 


2nd WING PENDING 



17.19 

16.80 


NACELLE VERTICAL 



TEST(HZ.) 19.14 

ANALYSIS(H2.) 19.08 


1ST WING CHORDWISE 



NACELLE ROLL 



29.20 


25.20 


29.93 


3rd WING BENDING 



TEST(H2.) 32.03 

ANALYSIS(H2.) 32.41 


1st WING TORSION 



46.39 

46.63 


LEGEND: 

NODE LINES 

TEST 

ANALYSIS 


FIGURE MEASURED A CALCULATED FPEOUENCIES « NODE LINES FOR 
WING (FULL) - NACELLE(NCMINAL) 


114 



1st WING BENDING 



ANALYSIS(HZ.) 5.59 


2nd V/ING BENDING 



14,45 

U.G4 


NACELLE SIDE BENDING 



15.62 

15.75 


NACELLE VERTICAL 



ANALYSIS(HZ.) 18.65 


1st WING CHORDWISE 



21.48 

22.87 


3rd WING BENDING 



27,73 


NACELLE ROLL 



ANALYSIS(HZ.) 30.26 


1st WING TORSION 



41.89 


LEGEND: 
NODE LINES 


TEST 

analysis 


FIGURE C 6 MEASURED S CALCULATED FPECUENCIES S NODE LINES FOR 
WING (FULL) - NACELLE! NOMINAL l-WINGLET (20 DEG) 
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APPENDIX D 


PR'CEDURE FOR MODIFYING STIFFNESS MATRIX 



^ k' 


APPENDIX D - Procedure for Modifying Stiffness Matrix 
The procedure is based on following assumptions: 


ORIGINAL P, .;. 

OF POOR QuAU: V 


1. The analytical mode shapes exactly match the test mode shapes. 

2. The frequecies for modes m + 1 through n are exact where m lowest 
frequencies are available from the test and n is the total 
degrees-of-freedom of the analytical model. 

3. The analytical mass distribution accurately describes the model. 






The pcocedure v/as applied >:o collect fhe stiffness matrix based on the clean 
wing vibration test. 





APPENDIX E 


SUMMARY OF EXPERIMENTAL RESULTS - HIGH SPEED MODEL 
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,U»WfiY OF EXPERIMENTAL RESULTS - HIGH SPEED MODEL 
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